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Abstract. We constraint on computer the best linear unbiased generalized 
statistics of random field for the best linear unbiased generalized statistics of 
an unknown constant mean of random field and derive the numerical gener- 
alized least-squares estimator of an unknown constant mean of random field. 
We derive the third constraint of spatial statistics and show that the classic 
generalized least-squares estimator of an unknown constant mean of the field 
is only an asymptotic disjunction of the numerical one. 



1. The best linear unbiased generalized statistics 
Remark. To simplify notation we use Einstein summation convention then 
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is given matrix. 

Let us consider the random field Vj] j £ Ni with an unknown constant mean 
m and variance a 2 its estimation statistics Vj and the variance of the difference 
Rj =Vj—Vj, where E{Vj} — E{Vj} = m, as covariance 

DHVj-Vj} = CoviM-VjXVj-Vj)} 

= Cov{Vj Vj } - Cov{V 3 Vj } - Cov{Vj Vj } + Cov{V Vj } 
Cov{VjVj} - 2Cov{VjVj} + Cov{VjVj} 
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and the linear estimation statistics (weighted variable) Vj = Y^i=i = 3 c 

i = 1 , . . . , n at j > n + 1 then 

D 2 {Rj} = Cov{V j V j }-2Cw{V j V j } + Cov{V j V j } 

= Var{V 3 } - 2Cov{J2^V t Vj} + G^E^OE^)} 

i i i 

= a 2 - 2j2^pov{ViVj} +J2J2 u) H Cov ( V ^ Vl } 
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(1) 



a z - 2a 2 \uj)p r] \ + (J 2 \u)p a uj) 
a 1 ± 2<T 2 w i j p i j =p a 2 Lu l J p u uJ l j , 



where pij] i = 1, . . . ,n is given vector of correlations and pa; i, I = 1, . . . , n is given 
(symmetric) matrix of correlations (see Appendix |^J) . 

The unbiasedness constraint (the first constraint on the estimation statistics) 

n 

E{R,} = E{Vj - Vj} = E{Vj} - E{Vj] = E{Vj} - E{w)Vi} = m - m]T u>) = 

i=i 

equal to 

n 

(2) £«} = / aiW } = a;j/„=l 

gives the first equation 
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The minimization constraint (the second constraint on the estimation statistics 
the statistics is the best) 



(3) 

where (JTJ) 



dD 2 {R 3 ) 



dujj 



±2a 2 p l0 =F 2a 2 p a J, J =F 2<j 2 f ll p 1 J = , 



D 2 {R J } = a 2 ± 2a 2 Lj} Pij T <? 2 u)paUj T 2a 2 (wj/« - l) M. 



produces n equations in n + 1 unknowns the kriging weights uij and a Lagrange 
parameter p* 
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this system of equations if multiplied by u)j 

UjPilUJj + Ujfil (ij = OJjPij , 

1 

and substituted into 

D 2 {R } = EWi-Vin- fW-Vi} 

£{[(V$-ro)-(V$,-m)] 2 } 
= E{[V 3 - m] 2 } - 2(E{V j V j } - m 2 ) + E{[Vj - m] 2 } 
= cr 2 - 2o- 2 \u)pij | + a 2 \uj l j p l iOJ l J | 

= cr 2 ± 2a 2 u>)p lj =F o 2 uj)pu<J j 

since variance of the (estimation) statistics is minimized 

Em-m} 2 } = CowKwjViXwjVi)} 

i I 
= <T 2 \u) Pl iJ 3 \ 

= T^Oj'jPilLOj 

(4) = T^{io} Pij - pj) 
gives 

(5) D 2 {R 3 } = £{[V- - V,} 2 } = E{[(Vj - m) - (Vj - m)] 2 } = a 2 (l ± (wjpy 
the constraints ([2]) and ([3]) produce n + 1 equations in n + 1 unknowns 
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2. The CLASSIC BEST LINEAR UNBIASED GENERALIZED STATISTICS of an 
UNKNOWN CONSTANT MEAN OF THE FIELD 

Remark. When we consider an independent set of the random variables Vf, i = 
1, . . . , ri with an unknown constant mean m and variance a 2 the best linear unbiased 
ordinary (estimation) statistics Vj = wjV* of the field Vj; jCt = l,...,n has the 
asymptotic property 

(7) lim £7{[u4Vj - m] 2 } = 

whilst for spatial dependence between random variables (the best linear unbiased 
generalized statistics) we get (see Appendix |Bj 

(8) lim lim E{[uAVi - m] 2 } = . 

n— > oo j—tao 

Due to different asymptotic limits between (J7J and ((8} the ordinary least-squares 
estimator of an unknown constant mean m of the field, the best linear unbiased es- 
timator of an unknown constant mean m of the field, can not be so easy generalized 
(like it was in past). 
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Let us constraint the best linear unbiased generalized (estimation) statistics Vj = 
UjVi of the random field Vj] j C i = 1, . . . , n, when for finite n and j — > oo the 
vector of correlations simplifies to 



£ -> 0" (or £ ^ 0+) 
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lim WjPij = Sulfa = £ 



then from |2j) 
(10) 

it holds flU 

(11) lim - ^] 2 } = lim ^(1 ± (wjptf + M j )) = ^(1 ± (£ + /zj )) 

for the co-ordinate independent statistics of an unknown constant mean of the field 
Vj with the constraint on (fTT]) 



lim E{[Vs - wjV^] 2 } = a 2 = E{[Vj - m] 2 } 



(12) 

given by constrained from (jlip 

(13) pj = -£ 

and from ([9]) the system of equations ([6]) 
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equivalent to 
and 
where 

w = 

nxl 

with the solution 
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of the classic best linear unbiased generalized statistics for finite n and j — > oo of 
an unknown constant mean of the field 

F'A^V 

(16) " 



lim w'V = 
j^oo F'A-^F 



where 



V = 



Vi 



v„ 



with constrained minimized variance of the best linear unbiased generalized (esti- 
mation) statistics ((4]) as its variance ffrom tflO]) and (fT3)) ) 

lim E{[(AVi - to] 2 } = lim T^Upij - /ij) = T<x 2 (£ - rf) = 



J-fOO 

then (fromflTH)) 



lim E{[w'V - to] 2 } = Tc 2 2£ 



j^oo F'A~ 1 F ' 

with the classic generalized least-squares estimator for finite n and j — > oo of an 
unknown constant mean to of the field 



(17) 



lim w v 

J-Voc 



F'A" 



F'A-^F 

based on observation v seen as outcome of V. 



3. The numerical best linear unbiased generalized statistics of an 
unknown constant mean of the field 

To remove the asymptotic limit of the classic best linear unbiased generalized sta- 
tistics for finite n and j — > oo of an unknown constant mean to = E{Vj} of the 
field Vj with the constraint (fT2|) 

lim E{[Vj - uWi] 2 } = a 2 = E{[Vj - to] 2 } , 

j->oo 

the best linear unbiased generalized (estimation) statistic of the field Vj ; j C i = 
1, 



, . . . , n at finite j > n 
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the negative correlation function with the parameter t — 182 + 1, . . . , 182 + 139 
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p(A y ) 
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for Aij = \i-j\>0, 
for = \i-j\=Q, 
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Figure 1. Variance of the numerical best linear unbiased gen- 
eralized statistics for finite n at finite j > n + 1 = f 82 + 1 of an 
unknown constant mean m — E{Vj} of the field Vj in units of 
the variance a 2 of the field computed 139 times for the negative 
correlation function (fT8"f with the parameter t = 182 + 1, . . . , 182 + 
139. 



was constrained (from ([5|)) on computer (139 times) for the numerical best linear 
unbiased generalized statistics for finite n at finite j of an unknown constant mean 
m = E{Vj} of the field Vj with the third constraint of spatial statistics 

(19) E{[Vj - c^] 2 } = a 2 = E{[V 3 - m] 2 } 
equivalent to 

(20) u)p VJ + M j = 

with constrained minimized variance of the best linear unbiased generalized (esti- 
mation) statistics (0]) as its variance (see Fig. [1} . 

Our aim was to derive for the negative correlation function (|18[) with the parame- 
ter t — 182 + 1 , . . . , 182 + 139 the numerical generalized least-squares estimator uj^Vi 
of an unknown constant mean m = E{Vj} of the field Vj in fact the proper best 
linear unbiased (generalized) estimator of an unknown constant mean m — E{Vj} 
of the field Vj given at finite j > n + 1 = 182 + 1 by numerical approximation to 
root of the equation (|20|) . This (co-ordinate dependent) generalized least-squares 
estimator wjfi was compared to the (co-ordinate independent) classic generalized 
least-squares estimator lim^oo u>)vi of an unknown constant mean of the field p7|) 

F'A -1 v 
lim w v = ; — ; — 

j^oo F'A _1 F 

based on the same observation an initial amplification vt = Vi , . . . , vis2 of long-lived 
asymmetric index profile recorded by 600 close quotes of Xetra Dax Index shown 
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Figure 2. Long-lived asymmetric index profile, Xetra Dax Index 
from 23 X 1997 up to 10 III 2000 (600 close quotes) the numerical 
generalized least-squares estimator uj*Vi of an unknown constant 
mean to = E{Vj} of the field Vj] j C i = 1, . . ., 182 (black dots) 
based on = V\ , . . . , «ig2 is compared for the negative correlation 
function (|18p with the parameter t = 182 + 1, . . . , 182 + 139 at finite 
j > n+1 = 182+1 to the classic generalized least-squares estimator 
linij_ ! . 00 ^JjVi of an unknown constant mean to = E{Vj} of the field 
Vj (grey line) with the same correlation function and based on the 
same sample. The classic estimator is the first approximation of 
the numerical estimator at final j = 577 for final t — 182 + 139. 
The dashed vertical line represents j = n = 182. 



in Fig. [2] then 



"182 



with the same correlation function (jl8|) . 

Since the classic best linear unbiased generalized statistics for finite n and j — >• co 
of an unknown constant mean m = E{Vj} of the field Vj with the constraint 

lim E{[Vj - com] 2 } =(T 2 = E{[Vj - to] 2 } , 

is an asymptotic disjunction for j — y co of the numerical best linear unbiased 
generalized statistics for finite n at finite j of an unknown constant mean m — 
E{Vj} of the field Vj with the constraint 



E{[V 3 ~Lo)V] 2 } = a 2 = E{[V 3 -m] 2 } 



LOjV l 



of an un- 
co of the 



then the correct classic generalized least-squares estimator linx,- 
known constant mean to of the field is an asymptotic disjunction for j 
numerical generalized least-squares estimator of an unknowm constant mean 
to of the field (see Fig. [5]). 
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4. Summary 

It was shown that the (estimation) statistics of the field Vj; j C i = l,...,n with 
an unknown constant mean to and variance a 2 

n 
i=l 

that assumes - the unbiasedness constraint ([2]) 

E{V j }-E{u i j Vi} = 
that assumes - the minimization constraint ([3]) 

dw) 

given by the kriging system of equations © 
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is the best linear unbiased generalized (estimation) statistics of random field Vj 
with minimized variance of the statistics 

(21) E{[u)Vi - to] 2 } = +V (u}/hj - Hi) 
and (minimized) 

(22) E{\yj - u^Y) = a 2 (1 ± + M j )) 
with the asymptotic property (Appendix [B|) 

lim lim E{[u)Wi - to] 2 } = 

and 

lim lim £J{[K--^yi] 2 } =a 2 

constrained once again from (|22[) on computer - the third constraint of spatial 
statistics 

E{[Vj - ojtyf] = a 2 = E{[V 3 - to] 2 } 

is the numerical best linear unbiased generalized statistics for finite n at finite 
j of an unknown constant mean to = E{Vj} of the field Vj with the numerical 
generalized least-squares estimator wji>j of an unknown constant mean of the field 
and its asymptotic disjunction for j — > oo the classic generalized least-squares 
estimator lim^oo uj^Vi of an unknown constant mean of the field. 
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Appendix A. The sign of the terms 

If for correlation matrix pu\ i,l = l,...,n that consists of unit diagonal ele- 
ments (see ([IE])) and non-positive off-diagonal elements holds 

oj)Puoj] < 

like at j > n+1 for vector p^ ; i = 1, . . . ,n that consists of non-positive correlations 
holds 

oj)p lj < 

then (P) 

D 2 {R J } = a 2 + 2a 2 u) PlJ - o 2 u)p u w\ 
for non-negative correlation function 

D 2 {Rj} = a 2 - 2a 2 uj) Pl] + a 2 uj)p u ^ 3 

for white noise 

D 2 {R j }=a 2 + a 2 u] Pa J j , 
where pu is the identity matrix. 

Appendix B. The asymptotic property of the best linear unbiased 

GENERALIZED STATISTICS OF RANDOM FIELD 

From the minimization constraint ([3j) 

PisUj + fu p.] = Pij 



P H pis = &\ , 

fll^J = 1 
n- =(/i//«)- 1 (/Wx-i) 



we get 
where 

substituted into ^ 
gives 
for j ->• oo ([9]) 

Pij = Lhi 

then the Lagrange parameter simplifies to 

lim pj = t-(fxiP H fu)- 1 

from the unbiasedness condition ([2]) it also holds 

lim v'pij = £,uj)f ll = £ 

then the minimized variance of the estimation statistics (j4]) 
E{[uty - m] 2 } = T° 2 {u)p % i - p]) 

simplifies to 

lim E{[u>iVi - m} 2 } = T <J 2 (t - Z + 

and © 

- W ^] 2 } = a 2 (l ± (cj} Pij + p> )) 
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simplifies to 

lim E{[Vj - ^Vi] 2 } = a 2 (l ± (£ + £ - (fii/Ffu J" 1 )) 

since 

lim (fulfil)' 1 = 

n— >oo 

and 

we get the asymptotic property of the best linear unbiased generalized statistics of 
random field 

lim lim E{[w)Vi -to] 2 } = 

n— >oo j— >oo 

and 

lim lim E{[Vj - uWi] 2 } = a 2 . 

n—tcx) j—>oo 
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